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ABSTRACT 

The spectral shape of the prompt emissions of gamma-ray bursts (GRBs) is typi¬ 
cally expressed by the Band function: smooth joining of two power-law functions for 
high- and low-energy regions. To reveal the origin of the Band function, we revisit 
the stochastic acceleration model, in which electrons are accelerated via scattering 
with turbulent waves in the GRB outflow. The balance between the acceleration and 
synchrotron cooling yields a narrow energy-distribution similar to the Maxwellian 
distribution. The synchrotron spectrum becomes consistent with the observed hard 
photon index for the low-energy region. On the other hand, the narrow electron en¬ 
ergy distribution contradicts the power-law spectrum for the high-energy region. We 
consider an evolution of the electron energy distribution to solve this problem. The 
turbulence and magnetic field induced by a certain hydrodynamical instability grad¬ 
ually decay. According to this evolution, the typical synchrotron photon energy also 
decreases with time. The time-integrated spectrum forms the power-law shape for the 
high-energy region. We discuss the required evolutions of the turbulence and mag¬ 
netic field to produce a typical Band function. Although the decay of the turbulence 
is highly uncertain, recent numerical simulations for decaying turbulence seem com¬ 
paratively positive for the stochastic acceleration model. Another condition required 
to reconcile observations is a much shorter duration of the stochastic acceleration than 
the dynamical time-scale. 

Key words: acceleration of particles - radiation mechanisms: non-thermal - turbu¬ 
lence - gamma-ray burst: general 


1 INTRODUCTION 

Gamma-ray bursts (GRBs) release isotropically equivalent energy of Eieo = 10®^-10®^ erg as gamma-rays within a time-scale 
of 0.1-100 s. The gamma-rays are considered to be emitted from relativistic jets with bulk Lorentz factor F > 100. The photon 
spectra are expressed with the conventional Band function (Band et al. 1993) with the parameters of the spectral peak energy 
Ep (typically 100 keV-1 MeV), low-energy photon index a (~ —1), and high-energy photon index fi (roughly between —2 
and —3). The short variability (1 ms-1 s) of GRB pulses is attributed to intermittent ejections of multiple emission regions 
from the central engine. In the classical internal shock model, this emission regions are shocked shells, where electrons are 
accelerated via the first-order Fermi acceleration (Fermi-I). However, the observationally typical low-energy index contradicts 
this scenario. Synchrotron emission from electrons injected via the Fermi-I process yields a = —1.5, because the strong 
synchrotron cooling leads to an electron energy distribution of N{Ee) oc E~^. 

Alternative ideas to resolve the low-energy photon index problem are photosphere models (Meszaros & Rees 2000; 
Giannios 2006; Pe’er et al. 2006; Asano & Meszaros 2013, and references therein), the Klein-Nishina effect on synchrotron 
self-Compton (SSC) process (Derishev et al. 2001; Bosnjak et al. 2009; Nakar et al. 2009; Wang et al. 2009; Bosnjak & Daigne 
2014, and references therein), sudden damping of magnetic fields (Pe’er & Zhang 2006; Zhao et al. 2014), and inverse Comp¬ 
ton models (Ghisellini & Celotti 1999; Stern & Poutanen 2004; Vurm & Poutanen 2009). In this paper, we revisit the second 
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order Fermi acceleration (Fermi-II) model of Asano & Terasawa (2009) (see also Mnrase et al. 2012). The Fermi-II process 
continuously accelerates electrons during the photon emission period. In this process, the resultant electron spectral index 
becomes harder than —2 (e.g. Schlickeiser 1984; Park & Petrosian 1995; Becker, Le & Dermer 2006; Stawarz & Petrosian 
2008; Lefa, Rieger & Aharonian 2011). As shown in Asano & Terasawa (2009), the acceleration time-scale required to balance 
with the synchrotron cooling is much longer than the shock acceleration time-scale with the Bohm limit; in other words, the 
electron mean free path in the Fermi-II models is much longer than the Larmor radius ru. The scattering efficiency required 
in the Fermi-II models is not so high compared to the Fermi-I acceleration. Asano & Terasawa (2009) showed that the Fermi- 
II model naturally reproduces the low-energy photon index of the GRB spectrum. The required turbulence to accelerate 
electrons may be induced via the Kelvin-Helmholtz instability in the shear flow (e.g. Zhang, Woosley & MacFadyen 2003; 
Mizuta & Aloy 2009) or at the boundary between the jet and cocoon (Meszaros & Rees 2001; Ramirez-Ruiz, Celotti & Rees 
2002). As radial modes, the Rayleigh-Taylor and Richtmyer-Meshkov instabilities are possible candidates to induce turbulence 
(Matsumoto & Masada 2013). In internal shock simulations with density fluctuation, the Richtmyer-Meshkov instability is 
excited (Inoue et al. 2011). 

Interestingly, recent studies with the Fermi-II process also agrees with broadband blazar spectra (Asano et al. 2014; 
Diltz & Bottcher 2014; Kakuwa et al. 2015; Asano & Hayashida 2015). The models reasonably reproduce very hard spectra 
seen in lES 0229+200 and 3C 279, and curved spectra for Mrk 421 and Mrk 501. Such results encourage us considering the 
Fermi-II process in GRB jets as well. 

While numerical simulations of the Fermi-II were done in Asano & Terasawa (2009) and Murase et al. (2012), we ana¬ 
lytically discuss the stochastic acceleration and emission in this paper. Whereas a simple Fermi-II model can easily produce 
observationally typical value of a ~ —1, we mainly focus on the high-energy slope of the photon spectrum. The simplest 
Fermi-II model implies a Maxwellian-like narrow energy distribution for electrons, which contradicts the power-law photon 
spectrum above £p. To resolve this problem, we consider temporal evolutions of the electron energy distribution, which may 
be regulated by the evolutions of the turbulence, mean magnetic field, and electron injection process. The resultant evolution 
of the photon peak energy would produce a power-law envelope curve in the photon spectrum. Our purpose in this paper is to 
find ideal evolutions of the model parameters to reproduce the GRB photon spectrum. Although the lack of knowledge with 
the turbulence property and evolution in GRB jets prevents the final confirmation of the obtained results, our study will be 
an important step to probe the possibility of the Fermi-II process in GRB. 


2 STOCHASTIC ACCELERATION AND COOLING 


Gamma-rays are emitted from a collimated outflow relativistically moving towards observers. Hereafter, we discuss the electron 
acceleration and photon emission in the plasma rest frame. The values in the following equations are defined in the rest frame 
except for photon energy; e' in the rest frame and e = Fe' in the observer frame. We start from the same Fokker-Planck 
equation in Asano & Terasawa (2009) as 


ON _ 



( 1 ) 


for the energy distribution of ultrarelativistic electrons N{Ee) in a certain volume. The energy diffusion coefficient Dee is 
phenomenologically assumed to have a power-law form 


Dee = KqE^. 


( 2 ) 


For simplicity, we consider only synchrotron radiation for the cooling processes for electrons. The energy loss rate for an 
electron of Ee = 'ye'meC? is expressed as 

Ecooi(7e) = -cyTcylUB, (3) 

where ut is the Thomson cross section, and Ub = B^/87r is the energy density of the magnetic field. The cooling time 


^c(7e) 


GnmeC 
o-TB^7e ’ 


(4) 


for gamma-ray emitting electrons is much shorter than the dynamical time-scale (see e.g. Asano & Terasawa 2009). Therefore, 
when the synchrotron cooling balances with the stochastic acceleration, the electron energy distribution can be approximated 
by the steady solution for eq. (1) with a pile-up feature (Schlickeiser 1984, 1985; Aharonian, Atoyan & Nahapetian 1986; 
Stawarz & Petrosian 2008; Lefa, Rieger & Aharonian 2011) as 


iV(Ee) 


SNtot (E,V 



3-9 


(5) 
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EeN(E,) 



where Ntot is the total electron number, F is the gamma function, and the cut-off energy is written as 

, 2„3 \ 1 /( 3 - 9 ) 


/ 67r(3 - q)KomlD \ 


( 6 ) 


As shown in Figure 1, the spectral peak energy in ElN{Ee)-'p\ot is sensitive to the index q. 

Adopting the synchrotron function J^syn(x) = Xwhere is the modified Bessel function of the 

second kind (Rybicki and Lightman 1979), the photon production rate can be written as 




where the typical photon energy for electrons of Eg is 

Substituting eq. (5) into equation (7), 


/ / T-i \ STvheB 2 

^typ(-£'e) — ~oTII 


8meC 


NRe') = 


Nu 


ttF 


(f5f ) h^El 


/ / 7? \ /’oo 

dFleexp - j , d^Ks/sR). 


Defining x = Ee/Eo and e'q = 4yp(-Eo), 

NRe') = 


Nu 


ttF 


(f5f) 


[ dxexp{-x^ '^) [ d^K^/sR) 

J Jx~‘^e' Ie'q 




ttF 


(f5f) h^E^, 


F, I 


(7) 

( 8 ) 

(9) 

( 10 ) 

( 11 ) 


Then, the spectral shape of the photon production rate is expressed by the non-dimensional function Fq{e/eo), which is 
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x' Fq(x) 



0.001 0.01 0.1 1 10 100 


Figure 2. The non-dimensional functions Fq that provide the spectral shape of the photon production rate Nry. 


Table 1. Xq = Cp/eo, the non-dimensional function lq(0,0) is defined by eq. (24). 


Q 

1 

3/2 

5/3 

2 

Xq 

2.5 

5.0 

7.2 

25 

7,(0, 0) 

1.1 

3.0 

5.3 

35 


obtained by integrating the fnnction over x and In Figure 2, we show the functional shape of Fq for various q. The 
horizontal axis is normalized by so- As shown in Figure 2, the spectral peaks in e^iVly(£)-plot are larger than so (Fritz 1989; 
Zirakashvili & Aharonian 2007; Lefa, Rieger & Aharonian 2011). Here, we define this peak as £p = XqSo, and summarize 


values of Xq in Table 1. 


We can choose fiducial parameters to agree with £p 


f Ko 


0.5 MeV as 
2 


500 J Vl02 s-i 
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B 


500/ Vl02 GeV^/® s-i 


Ko 


_L^ _ 

500/ Vl02 GeV^/^ s-i 


4/3 


,104 G 

3/2. ^ 

Vl04 G 
B 


F 


Ko 


102 GeV s-i / V 104 G 


104 G 
/ B 


^ MeV, 

for g = 2, 

(12) 

^ MeV, 

for q = 5/3, 

(13) 

-5/3 

MeV, 

for q = 3/2, 

(14) 

^ MeV, 

for g = 1. 

(15) 


Note that the typical electron energy to emit 0.1-1 MeV photons is ~ GeV in the source rest frame for F = 500 and 
B — 10^ G. The acceleration time-scale face oc implied above 10“^ s for GeV electrons) is much longer than 

the gyration period 2'KEel{eBc) ~ 7 x 10“®(i?e/GeV)(i3/10^ G)“^ s. The Fermi-II model implies that the energy gain 
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per one scattering is AEe/Ee ^ where Peg is the effective turbulence velocity normalized by c. The energy fraction 
of the turbulence to the bulk jet energy is ~ /3|g. Though the turbulence velocity may be slower than the sound velocity 
in the relativistic plasma {p'^g < 1/3), the high radiative efficiency (e.g. Zhang et al. 2007) may require p^g ^ 0.1. This 
situation is similar to the standard shock acceleration models, in which a mildly relativistic shock is presumed. In both 
the cases, the energy source of the non-thermal electrons is the energy of mildly relativistic protons. Here, we assume that 
the turbulence has a characteristic scale, namely typical eddy size ledd- If the turbulence spectrum ranges in significantly 
short wavelengths comparable to the Larmor radius rg ~ 300(i5e/10® eV)(B/10'* G)“^, the mean free path is written as 
^mfp ~ (Blandford & Eichler 1987). Inside the eddies P^g <C 0.1, while only electrons moving across eddies are 

scattered via nonresonant mirror interaction (transit time damping) arose from compressible (acoustic) modes (Ptuskin 1988; 
Cho & Lazarian 2006; Yan & Lazarian 2008; Lynn et al. 2014) with P^g ~ 0.1. When the diffusion time ledd/GGfp) is much 
longer than the eddy time-scale tedd = ledd/{cPeg), the time-scale t^ccPeg may correspond to the eddy time-scale so that 
Zedd ~ Ctacc/3efi ~ lO’^(tacc/10“^ s)(,0efj/0.1 )^/^ Cm. Finally we obtain Dee ~ PlgEl/Ugg oc Plg/l^dd- 

However, there may be a transition scale at the wavelength where the kinetic energy is comparable to the magnetic 
energy. At the shorter wavelength than the transition scale, the power-spectrum of the turbulence is steeply damped. Then, 
the mean free path of electrons may be elongated as long as the eddy scale. In this case, we obtain Zedd ~ ct^ccPeg ~ 
3 X 10’’(tacc/10“^ s)(/3eff/0.1) cm and Dee oc Pig/Udd- 

The above qualitative discussion is still insufficient to conclude the acceleration property in the turbulence. In the following 
discussion we maintain the phenomenological parametrization for the energy diffusion coefficient with arbitrary q and Kq. 


3 DECAYING TURBULENCE 


If the parameters in equation (11), such as eo etc., remain constant during the emission period, the photon spectrum for 
observers is just a blue-shifted spectrum of Afy (e'). As shown in Figure 2, the spectral shape with a sharp cut-off may contradict 
the observed Band spectra especially for smaller q. However, the efficiency of the stochastic acceleration is expected to decay 
with time. Then, the typical photon energy eo would shift to lower energy. Such an evolution of the photon emission spectrum 
can produce a high-energy power-law part in the time-integrated spectrum. 

Since the time-integrated photon spectrum is written as 


NRe') 



(16) 

(17) 


the temporal evolutions of the two quantities, eo and Ntot/Eo, determine the final photon spectral shape. Let us assume that 
all the parameters evolve following power-laws of (t -|- ftbi) and the turbulence suddenly disappears at t = tf for simplicity. 
The duration t{ is defined as the energy transfer time-scale from the turbulence to electrons. The parameter ttbi, which 
may correspond to the production time-scale of the turbulence, is introduced to avoid divergence at t = 0. Hereafter, the 
combinations of the parameters are assumed to behave as Ntot/Eo oc (t -I- ttbi)“ and e'o oc BEq oc {t + ftbi)”** with indices a 
and b. Here we explore ideal pairs of the phenomenological indices a and b to reproduce typical GRB spectra. 

As will be discussed below, the index a may be < 2. In Figure 3, we plot the photon spectrum N-y{e) with q = 5/3, a — 2, 
and tf = lOOftbi for various b. The photon energy is normalized by £o,i, which is £o at f = 0. As b becomes large, the faster 
decay of eo makes a lower spectral peak energy. Since f dtJVj(eo) oc (t -I- ttbi)^^“, the photon index parameter in the Band 
spectrum would be approximated as 




1 A n 
6 


(18) 


Therefore, as shown in Figure 3, the parameter values of b between 1 and 1.5 are favorable for a = 2 to adjust the index 
between —2 and —3. This result is common irrespective of the value of q\ the spectral shapes for other q are similar to each 
other. Fspecially, the case of fe = 1 in Figure 3 is well fitted with the typical Band parameters {a = —1 and p = —2.5) below 
10£p. The obtained Band parameter £p is slightly larger than the final value of Xqeo ~ 7.2{t + ttbi)~^£o.i — 0.07£o,i, but the 
difference may be hard to be seen as shown in the figure. The spectral curvature around £p leads to the difference between 
the Band parameter P = —2.5 below 10£p and the asymptotic value —(1 + a)/b = —3. Note that, however, in most cases the 
photon statistics above 10£p (typically a few MeV) is insufficient to detect the spectral break or curvature. 

Here we attempt to associate the phenomenological parameters a and b with the evolutions of the turbulence and electron 
injection. Initially the diffusion coefficient and magnetic field are larger than the values assumed in eqs. (12) -(15). Then, 
the acceleration time-scale is gradually elongated, and settles down to the value suited for £p ~ 0.5 MeV at t = tf. Let us 
write Ko oc t~‘‘ and B oc t~™, Eq oc Then, e'o oc BEq oc i-™+2(2™-s)/(3-q)^ total number of electrons may 
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£^N^(£) 



0.001 0.01 0.1 1 10 100 1000 


s/eo,i 

Figure 3. Photon spectrum with q = 5/3, a = 2, and tf = lOOitbl for various b. The dashed line is the Band function with a = —1, 
/3 = —2.5, and Sp = 0.13eo,i. 


increase following the injection rate If we express this as Ntot oc f", tN^le'q) oc iNtotE^ ^ oc As a result, 

the parameters a and b are written with the indices defined above (s, w, and n) as 


a = n — 2 


2w — s 
3-q ’ 


b = w — 2 


2w — s 
3-q 


(19) 


Although the parameters s, w, and n are highly uncertain, some discussion of the interaction between the turbulence 
and electrons may provide insight into the required condition, (1 -\- a)/b = 2-3, in our Fermi-II model. The average energy 
gain of electrons per one scattering by waves is proportional to ~ max /3 a), where /3w and Pa are the turbulence 
velocity and Alfven velocity normalized by c, respectively. In our case, the required electron energy density is higher than 
the magnetic field density so that /3|v > /3 a is likely. In the simulations of Inoue et al. (2011), when a relativistic shock 
propagates inhomogeneous medium, the Richtmyer-Meshkov instability is excited and the turbulence starts to decay after a 
few eddy-turnover times. If we adopt those results, 0^ oc (7k oc and 5B^ oc Ub oc As we discussed in the last part 

of §2, the diffusion coefficient may behave as Dee oc /f^/ledd (short-mfp model) or /3^/(edd (long-mfp model). The long-mfp 
model with a constant (edd and the simulation result of Inoue et al. (2011) lead to g = 2, s = 1.3. Since the average field 
B oc 5B in this highly turbulent plasma, we obtain w = 0.35. The above values imply b = 1.55 and a = n -|- 1.2. Therefore, 
the constant injection rate n = 1 seems favorable, while the short-mfp model requires an ad hoc value of n ~ 3. 

Other recent MHD simulations (Zrake 2014; Brandenburg, Kahniashvili & Tevzadze 2015) show inverse transfers from 
small to large turbulence-scales in spite of the absence of magnetic helicity. In such cases, the evolutions of the turbulence scale 
and energy densities depend on their initial conditions (Olesen 1997). In the relativistic simulation in Zrake (2014), the longest 


1 


Note that the electron number was assumed to be constant in the simulations of Asano & Terasawa (2009) 
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wavelength evolves as fcivi oc and Ub oc while the non-relativistic simulation of Brandenburg, Kahniashvili & Tevzadze 

(2015) shows kM oc and Ub oc . The indices of the power-spectral density for k > Am are 2 in both the two simula¬ 
tions. The simplest analytical discussion in Son (1999), considering the evolution of patchy turbulence, concludes ku oc 
and Uk oc [/b oc which are not far from the simulation results. From the derived indices ui = 3/5 and s = 8/5 for 

the long-mfp model, finally we obtain the phenomenological indices a = n -|- 4/5 and b = 7/5. If the electron injection rate 
is constant (n = 1), the photon index becomes an ideal value p = —(1 a)/h = —2. This primitive model of the turbulence 
evolution also predicts a reasonable range of the photon index. The electron acceleration affects the decay of turbulence, 
which modifies the estimate of the photon index. More quantitative discussion based on numerical simulations is necessary to 
determine the value p. On the other hand, the short-mfp model requires n = 3 again. 

Although the above discussion does not always validate the Fermi-II model, we do not find any discouraging evidence so 
far. The next step is a study to verify the fundamental acceleration process in such turbulence, but this is beyond our scope 
in this paper. 

We have neglected the effect of inverse Compton scattering. However, most of GRBs do not show the extra component 
that can be interpreted as the inverse Compton component (Ackermann et al. 2012). Even when the magnetic field is weak 
compared to the electron energy density, the lack of the inverse Compton component can be explained by the Klein-Nishina 
effect (e.g. Asano & Meszaros 2011). The numerical simulations of Asano & Terasawa (2009) also shows that the inverse 
Compton emission does not greatly affect the synchrotron spectrum. 


4 SECONDARY ELECTRON-POSITRON PAIRS 


The continuous electron acceleration increases the gamma-ray energy emitted from one electron so that the required electron 
number is less than the standard shock acceleration models, in which all electrons are assumed to be accelerated. In our model 
there should be a condition for electrons to enter the acceleration process. The injection mechanism is uncertain similarly to 
that in the shock acceleration. The electron injection mechanism may be sub-shocks (Narayan & Kumar 2009) or magnetic 
reconnection in the turbulence. At least, the electron should be relativistic to be in the acceleration process we have assumed. 
Hereafter, we assume the threshold energy of electrons/positrons, 7 thJTieC^, to be accelerated by turbulence. 

We have assumed a gradual increase of the electron number by a stable electron injection. However, high energy photons 
produce secondary electron-positron pairs via 77 -absorption. Those pairs can be also accelerated by turbulence and emit high- 
energy photons. If such pairs dominate the lepton number, the non-linear growth of the electron/positron number occurs, 
which results in a sudden increase of pairs. Those non-linear evolution of the electron/positron number is not favorable to 
produce a Band-like spectrum. In this section, we estimate the contribution of the secondary pairs under the situation we 
have assumed. 

First, we obtain the required number of electrons in the Fermi-H model. In our model, following the power-law evolution 
of (1 -\-t/Uh\)~^, the typical photon energy decays from £o,i to eo,f. Denoting the values of Atot and Eq at f = U as Atot.f and 
Eo,f, respectively, the total photon energy in the comoving frame is written as 


TP' /* , / /^r / /\ v3e Trifle Atot, 
Eph = dee NRe ) = --^ 

f f 

A''"/ 


1 + 


tthl 


Ue'li 


rtf/hbi 


dxy{l -I- xRFg 


(1 + ®)' 


( 20 ) 


The peak energy for observers is written with the final value of eo as ep ~ FA^eo f. Then, we obtain 


p ' _ 


1 -f 


Uhl J 


Substituting Eo,f/(meC^) for 7 e in eq. (8), we obtain 


gQ.i ^ / uV i-KheBf f 1 y 

E^j \ tthl) SrUeC VrrieC^y 


( 21 ) 


( 22 ) 


as well, where Bf is the final value of B. Adopting those equations, the isotropically-equivalent total photon energy in one 
pulse is expressed as 


£^ph = F J dee'N^le) 


3^36® 


8 F 


(fe?) 


hnif 


— Ntot,iB[tf^Ig{a,b), 

r >2 


(23) 
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where the non-dimensional function 

/ \ ^^ r /^tbi / 

s (l + ^j Jivl + 

for 1 -f a — 26 > 0, and the final approximation has been obtained with tf/ttbi S> 1. The values of lg(0,0) for various q are 
tabulated in Table 1. 

While the observables are -Bph and £p, the model parameters to determine the total photon energy are Wot.f, S{, and 
tf. The time-scale t{ can be considered as the dissipation time-scale of the turbulence into electrons. Alternatively, as will be 
explained below, the injection of secondary electron-positron pairs may control the duration t{. The time-scale tf can be shorter 
than the dynamical time-scale tdyn = R/{cT) as supposed in Asano & Terasawa (2009). Note tdyn = Ttv = 50(r/500)(tv/0.1 s) 
s, where tv = R/{cr^) is the variability time-scale for observers. But tf may be longer than the cooling time-scale 


(24) 

(25) 


tc(£p) 


Ott 


The total photon energy is written as 


AeTmeC 


= 2.6 


—) 

500 J 


1/2 


Bf 

10-* G 


-3/2 


(o.5 Mev) 


- 1/2 


Sph - 1.9 X 10® 


Atot,f 

To^ 


f Bf 


tf 

0.1 s J VlO-* G 


( -) 

V0.5 MeV/ 


erg, 


(26) 


(27) 


for q = 5/3 with 7,(0,0). The required electron number Atot.f to achieve ~ 10®^ erg is much less than the number in the 
classical shock acceleration model by a factor of ~ tc/tf. The assumed time-scale above implies a significantly longer scale 
than the eddy scale estimated in §2 as ctf ~ 1004dd. 

Let us estimate the number of the secondary pairs. We approximate the photon-density spectrum for e' > ej, by a 
power-law 

P 


/ /\ '‘-0 / t 


(28) 


Although the result will be independent of the width of the emission region, here we assume the width as ctf/3. Then, the 
normalization in eq. (28) becomes 


_ -3(/3 -b 2)Epf, 

° AnR^cttep 

In this case, the optical depth for 77 -absorption is written as 


ctf ( e'e'p ^ 


~ O.lcTTno — 

(Asano & Takahara 2003, and supporting material in Abdo et al. 2009). This results in 

i 2x^1 “(i® + 2)7?ph „i+/3 / 7ep 

r,,(7mec):i.0.1aT T ^ 

When r.y-y( 7 thm-eC^) <C 1, the number of the secondary pairs above 'yth'rrieC^ is independent of 7 th as 


A‘±(7 > 7th) 


r.y^ (7thmeC^)47r7?^ — 7thmeC^n(y (7thmeC^) 


— 0 2)^CrT lHeC 7?pt, ^2(3-2 


47rc^tH 


TTleC^ 


-2/3-1 


For P = —2.5, the number becomes 
iV±(7 > 7th) 


2.8 X 10^ 


0.1 s 


E, 


1051 erg/ V500 


ph 


( " 1 - ) 

V0.5 MeV/ 


(29) 

(30) 

(31) 

(32) 

(33) 

(34) 


which is smaller than the number Atot.f adopted in eq. (27), but of considerable quantity. If the condition N± -C Atot.f is 
critical for the model, we need to adjust the model parameters, e.g., increase of Atot.f by a smaller Bftf or decrease of N± by 
a larger T. As shown in Figure 3, however, the photon spectrum may not be the extension of the single power-law as far as 
the energy of FqthmeC^ ~ 26(r/500)(7th/100) GeV. In such cases, the number of the pairs is smaller than the above estimate. 
The numerical simulations in Asano & Terasawa (2009) is a case where the contribution of the secondary pairs is negligible, 
while the pairs injected via the hadronic cascade dominates in the model of Murase et al. (2012). 

Another possibility is that the sudden increase of electrons/positrons due to the non-linear effect triggers the damping 
of the turbulence. As long as the total energy in the turbulence is finite, an explosive increase of the lepton number causes 
a strong dissipation of the turbulence. In that case, the average energy of the electrons/positrons in the acceleration process 
also suddenly decreases. This leads to inefhciency of gamma-ray emission. Such a non-linear effect between the turbulence 
and secondary pairs may finally determine the duration tf. 
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Figure 4. Observable photon spectra for various jet opening angles. 


In conclusion, the turbulence decay time-scale tf, which depends on the electron injection rate including the secondary 
pairs, should be much shorter than the dynamical time-scale in the stochastic acceleration model. Note that this time-scale 
is still much longer than the presumed energy-transfer time-scale from protons to electrons in the shock acceleration models. 


5 OFF-AXIS CONTRIBUTION 

In this section we switch a topic to focus on the low-energy portion of the Band spectrum. One may consider that the hard 
electron spectrum leads to the “synchrotron limit” a = —2/3 for the photon index. Such a hard spectral index is not frequently 
detected. However, the best model in Figure 3 (6 = 1) is well fitted with the typical Band parameter a = —1. The difference 
between the model and Band function is hard to be recognized observationally. Moreover, the contribution of the off-axis 
emission can soften the photon spectrum for observers. The dominant contribution to photons detected by an observer is 
emitted from the jet surface within the opening angle (relative to the line-of-sight) of l/F because of the relativistic beaming. 
However, the emission from 9 > 1/V can modify the spectrum as shown in Figure 4. 

Here, we have assumed that photons distributing isotropically in the rest frame are emitted from a sphere expanding with 
r = 300. The photon spectrum in the rest frame has the form of Fq with q = 5/3. We have numerically integrated photons 
emitted from the surface within 9 < 9j taking into account the relativistic beaming effects; photon energy e' = er(l —/3sh cos 9), 
solid angle dQ, — r^(l —/3sh cosphoton number dN = dN', surface dS — dS', angle cos9' = (cos0 —/Ish)/(1 —/3sh cos0), 
and dN'/{dtl'dS') oc | cos0'| (see Asano & Meszaros 2011), where /3sh = \/l — l/F^. 

The well-known photon index —2/3 for the synchrotron function (Rybicki and Lightman 1979) is obtained from Fq{x) 
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in the limit of x "C 1. Even for the range of 10“^ep“10”^ ®ep, however, the functions are fitted with a photon index between 
—0.84 and —0.79. This is seen in the low-energy slope for narrow jet opening angle of 6j = O.l/T or l/T in Figure 4. For the 
large opening angle of lO/F, the off-axis contribution slightly enhances the low-energy flux, which results in a photon index 
of ~ — 1. This slight modification obtained numerically may be hard to be treated analytically. If the jet has angle-dependent 
spectrum like the model of Lundman, Pe’er & Ryde (2013), various photon index would be obtained. 


6 SUMMARY 

We have discussed the stochastic acceleration by turbulence in GRB jets. Since the energy source is the turbulence, mildly 
relativistic turbulence is required to achieve a high radiative efficiency. When the electron spectral distribution is regulated 
by the balance between stochastic acceleration and synchrotron cooling, the temporal evolutions of the two combinations of 
the phenomenological parameters, Ntot/EQ and BEq, determine the final photon spectrum. If we choose an ideal parameter 
set, the photon spectrum can be fitted with the typical Band parameters a ~ —1 and P ^ —2.5. The off-axis contribution can 
soften the low-energy photon index further. The hard-to-soft evolution of photon spectrum due to the decay of the turbulence 
in our model may agree with the observed signatures (Norris et al. 1986; Bhat et al. 1994; Ford et al. 1995; Burgess et al. 
2014). The duration of the stochastic acceleration should be much shorter than the dynamical time-scale for hducial values 
of the magnetic field. The injection of the secondary electron-positron pairs may control the duration time-scale. 

The required parameter set is translated into the temporal indices for the evolutions of the diffusion coefficient, magnetic 
field, and total electron number. The excitation mechanism and evolution of the turbulence in GRB jets are highly uncertain. 
The verifications of the turbulence evolution and details of the acceleration mechanism are beyond the scope of this paper. 
However, the implication obtained from the recent numerical simulations of decaying MHD turbulence (Inoue et al. 2011; 
Zrake 2014; Brandenburg, Kahniashvili & Tevzadze 2015) seems encouraging for the required parameter evolutions. We can 
expect that future simulations of electron acceleration in turbulence, including the cases with magnetic reconnection (Hoshino 
2012; Kagan, Milosavljevic & Spitkovsky 2013; Dahlin, Drake & Swisdak 2014; Guo et al. 2014, 2015), reveal the excitation 
mechanism and initial condition of the turbulence to reproduce the GRB emissions. 
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